A mechanism of orbital reconstruction at the interfaces of transition metal oxides 
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Orbital reconstruction at interfaces between YBa2Cu306 and SrO-terminated SrTiOs is studied 
using local spin density approximation (LSDA) with intra-atomic Coulomb repulsion (LSDA+(7). 
The change of population of interfacial Cu 3d orbitals results in stabilization of a new oxidation 
state 3d* which involves an additional modification of orbital occupancies in the nearest SrO and 
Ti02 layers. We find that an increase of electron charge in Cu 3d^2_y2 states counterbalances a 
depopulation of 3d3^2_^2 orbitals which induces, on account of the onsite Coulomb repulsion U, a 
splitting of 3d^^2_^2 states at Cu02-SrO interfaces. 

PACS numbers: 74.78.Fk,74.78.-w,73.20.-r 
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I. INTRODUCTION 

Interfaces between transition metal oxides attract 
much attention due to their technological and fundamen- 
tal importance. In a complex transition metal oxide sys- 
tem, new electronic states at the interfaces between dif- 
ferent components can bring dramatic changes into its 
physical properties. The prominent examples are metal- 
lic mixed- valence heterostructures of SrTiOa and LaTiOs 
and the quasi-two-dimensional electron gas at the inter- 
face between two insulators SrTiOa and LaAlOa with 
recently discussed superconducting propertiesi^'^. 

In the transition metal oxide heterostructures, the 
new interfacial electronic states can be driven by the 
mechanisms of electronic reconstructions^. These mech- 
anisms involve charge compensation of interface polarity 
or charge self-doping due to interface chemical abrupt- 
ness which is one of the ways to maintain the overall 
electrostatical neutrality of the system^,. Besides the well 
understood scenarios of interface charge doping, other 
mechanisms of interface reconstructions can play a deci- 
sive role in the unusual electronic properties of transition 
metal systems. Recently the so called orbital reconstruc- 
tion has been shown to occur at the interfaces between 
cuprate and manganite films which results in changes of 
3d orbital occupancies of the interfacial copper ions^. 

The first evidence for the interface orbital reconstruc- 
tion came from resonant x-ray spectroscopy studies of 
YBa2Cu307/LaCaMn03-heterostructures^. In these ex- 
periments, the x-ray spectral intensity at the interface 
between YBa2Cu307 and Lao.67Cao.33Mn03 is found to 
be similar for the photon polarization in (x, ?/)-plane and 
in z direction, which is interpreted in terms of the appear- 
ance of hole charge in 3d3j,2_^2 orbitals. This is in strong 
contrast to the standard hole doping mechanism where 
the holes occupy O 2p orbitals of the Zhang-Rice singlets 
in Cu+^Oj^^ planes of YBa2Cu307_<5. The 'id^.z^^r'^ or- 
bital reconstruction has been explained by a formation 
of a new "extended molecular orbital" at the oxide in- 
terfaces in which a substantial amount of hole charge is 
distributed over the Mn, Cu and O orbitals. As is shown 
in Refs. dH, the high Mn charge -t-3 results in an electron 



transfer from Mn to Zdx2_y2 orbitals of Cu, whereas the 
Cu 3d3z2_r2 states become partially occupied. Moreover, 
in such an electron doped system, the coexistence of two 
partially occupied eg orbitals leads to the formation of 
Zhang-Rice triplets centered on Cu with the total spin 
S ~ \. In contrast, for the hole-doped interfaces, the au- 
thor in Ref. 8 expects a rather standard — y^-singlet 
character of the electron configuration with a stabiliza- 
tion of triplet states when the hole density is increased. 
As is shown in Ref. d, the singlet character of the electron 
state can be possibly changed due to structural Jahn- 
Teller distortions. 

In the present work, we explore theoretically a mecha- 
nism of the orbital reconstruction occurring through the 
splitting of 3(i322_r2 I and Sd^^^^r^ i^ bands by the inter- 
face charge redistribution. As an example of a system 
where the interface-caused splitting of 3d3z2_^2 can be 
obtained within the density functional theory, we con- 
sider a supercell with a Cu02 layer directly deposited on 
the top of a SrTi03 slab terminated by SrO, a model 
oxide interface recently studied in Refs. lolfiol and shown 
schematically in Fig.[TJ In this model system, the electro- 
static polarity of the interfaces between Cu02 and SrO 
could be compensated via the interface self-doping by 
holes. Like in the bulk YBa2Cu307_5, in heterostruc- 
tures the interface self-doped holes in the Cu02 layers 
are accumulated predominantly in the O 2p orbitals. Al- 
though the charge self-doping mechanism appears to be 
a common scenario in a wide variety of polar oxide inter- 
faces, a close examination of the electronic structure in 
the Cu02/SrO region unexpectedly shows that a mech- 
anism different from the electrostatic self-doping plays a 
key rolei^. 



II. INTERFACE CRYSTAL FIELDS AND 
STRUCTURAL RELAXATION 

As appears from the earlier LDA+J7 studies^' the 
relaxation of the interfacial distance h between Cu02 
and SrO planes results in the optimized value h = ho = 
1.83 A which is by 0.5 A smaller than the distance be- 
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FIG. 1: Scheme of a polar Cu02/SrTi03-superlattice where a 
STO-layer is terminated by a SrO-plane. The right side shows 
a structural configuration which appears at the interface. 



tween the apical oxygen 0(2) and Cu in Cu02 planes 
of the bulk YBa2Cu3 06. For larger interface distances 
h > 2.2A , the compensation of the polarity proceeds 
through a standard self-doping mechanism of the elec- 
tronic reconstruction and leads to an extremely high self- 
doping level nji = +2 holes/per unit cell which enforces 
charge neutrality (cf. the formal valency of each layer 
in Fig. [U to the left of the block scheme). As can be 
seen in Fig. [21 the holes are located in 0(1) 2px,y or- 
bitals (top panel, dashed arrows), in close analogy to the 
Zhang-Rice singlets described by a 3cP configuration in 
the Cu02 planes. 

In contrast to the unrelaxed interfaces, for the opti- 
mized distances h ^ the new (orbital) mechanism 
of the compensation of the interface polarity becomes 
equally important. This new mechanism is associated 
with a change of the occupancy of one of the ?id-^z2_j.2 
orbitals from one to zero due to a "displacement" of this 
band above the Fermi level (see Fig. [3]). From the point 
of view of classical chemical valence theory, the resulting 
electronic configuration can be formally described as a 
"3(i*" state. In Fig. [21 the charge and orbital reconstruc- 
tions caused by the -f2e-doping level are schematically 
shown for comparison. Formally, a Cu Sd^ configuration 
corresponds to a change of the valency of Cu from +2 
to -|-3 and is typically considered as unstable in the lit- 
erature. The expectation of such an instability was the 
reason for a preference for the standard self-doping re- 
construction mechanism in our earlier discussion of the 
Cu02/SrTi03 system in Ref. 9. However, a further close 
and careful examination of the electronic reconstruction 
in this system shows a more complex scenario which is 
beyond a simple integer-like change of the electron charge 
on the copper ion (Fig. [21 bottom panel). In this mech- 
anism, the modification of the oxidation state of Cu can 
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FIG. 2: Schematic view of the electron configurations which 
result from hole doping. The top and bottom right panels 
show the electronic and orbital reconstructed configurations 
with black arrows denoting the electronically fully occupied 
(bold arrows) and partially occupied (dashed arrows) orbitals. 
The configurations "3ci*" and "3d®" formally describe the 
states 3d* and 3d® with partially occupied Cg orbitals. 



be described by two steps, namely (i) by a depletion of 
electron charge from a ?id^^2_j,2 band and from the 2pz 
band of 0(2) (SrO), and (ii) by an increase of electron 
density in the idx^-y'^ orbitals hybridized with 0(1) 2px^y 
states. These two steps provide jointly a decrease of the 
net electron charge on Cu which is an additional way to 
compensate the interface polarity. Consequently, in the 
relaxed structure the electrostatic neutrality is achieved 
by a combination (a) of the standard self-doping mech- 
anism with the self-doped hole charge n?, = 1 on O 2p 
orbitals (Fig. [21 top panel), and (b) of the 'id^z2_j.2 or- 
bital reconstruction which leads to a decrease of the local 
electron charge on Cu (Fig. [21 bottom panel) . 

Let us estimate the crystal field which would be suffi- 
cient for stabilization of the new oxidation state formally 
denoted here as extended Cu 3d^. Using the notations 
{t, i} for the majority and minority spin states, the 3ci® 
state can be described as ?>d^^'^^ configuration where the 
occupied states Zdx2_y2 and 3d322_^2 ^ have a character 
of a Zhang-Rice triplet, similar to the results of Ref. d. 
In the relaxed structure, this configuration corresponds 
to the minimum of the total energy. In the LDA-I-J7 cal- 
culations, the energy gain AEtot = Etot{hi) - Etot{ho) 
due to the relaxation of the interface distance from 
h = hi — 2.2A to h — ho approaches 1 eV. It is re- 
markable that in the extended Cu 3d^ configuration, the 
main contribution to AEtot is provided by the reconstruc- 
tion of the electronic states in the nearest Cu02 and SrO 
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TABLE I: Band centers of gravity (in eV) for different values 
of interfacial distance h. Here, £2,2 refers to the 2pz orbital 
of 0(2). 



h (A), a 



1.83 (T) 
1.83 (i) 
2.20 (T) 
2.20 (i) 



-5.03 
-1.07 
-4.68 
0.97 



-5.23 
-4.84 
-3.60 
-3.16 



£2,z 



-2.98 
-2.62 
-1.93 
-1.82 



FIG. 3: Partial Cu 3d orbital density of states at the inter- 
face between CUO2 and SrTiOs terminated by SrO (LSDA-|-f/ 
studies) . The zero of energy is at the Fermi level. The left and 
right panels correspond to the case of unrelaxed {h = 2.2 A) 
and optimized (h = 1.83 A) interfacial distances, respectively. 

planes 

AEtot = Ae322_^2 + Ae,j.i_yi 

+Ae^,+y, + Ap^, (1) 

where A£3j,2_^2, A£a;2_j,2, and Asxz+yz are the changes 
of the Sd orbital energies; Ap^ refers to the shifts of the 
energies of the 2p bands of oxygens 0(2). In our calcula- 
tions, the band energies are represented by their gravity 
centersiiii^. Here a gravity center Eg of a band with a 
density of states p{E) is calculated from the expression 
Eg = J dEEp{E)/J dEp{E). As the displacements of 
the gravity centers of the 2p states of 0(1), 2px and 2py 
states of 0(2) and 3dxy bands of Cu are negligibly small 
as compared to the splitting of eg, 3dxz and 3dyz bands, 
they are not included into ([1]). 

In the SIC LBA+U approachi^, the Sd^z^^r^ orbital 
energies in the optimized (z/ = 0) and unrelaxed {v — 1) 
system can be represented as 




Here Ueff = U — J and a is the spin index. In our 
calculations, wc consider the local Coulomb corrections 
for Cu 3d orbitals U = 8 eV, and the Hund coupling 
J = 0.8 eV. In the unrelaxed structure, the states d^z^-r^ 
are filled by electrons: n}^^2_^2 g. ~ 1- In contrast, in 
the optimized system, the orbital reconstruction leads to 

"sz^-rM ~ ^ ^'^'^ "-Sz^-r^.i ~ 0- The quantities £322^^2 
are the corresponding LDA orbital potentials which are 
assumed to be independent of the spin direction. As 
a consequence, the total id^z^-r^ energy shift for the 
relaxed structure is given by 

A£3z2_r2 = Ae3^2_r2 jn322_r2,| -f £3^2 _^2 | ^3^2 _r2 ^ ^ 

where the potential difference A£i°2 „2 — £^^2'' 2 — 



£^\2_j,2 is determined by the change of the orbital crystal 
field due to the interfacial relaxation. 

From (21) , the change of the splitting of the spin bands 
£322_^2 I and £322_^2_| due to structural relaxation can 
be described as 

A£322_r2 - A£322_r2 

= C^e//(nL2-r2,i - "L^-r^a) ~ Ueff- (4) 

where n3^2_^2 ^ is the occupation of the (3z^ — r^, |) 
orbital in the optimized {v — 0) and unrelaxed {v = 
1) system. The last result demonstrates that the on- 
site Coulomb repulsion controls the spin band splitting 
and the consequent enhancement of the local magnetic 
moments. 

As the width of the id^z'^-r'^ bands exceeds the shift 
A£322_r2 caused by the orbital reconstruction, the use 
of the centers of gravity to estimate the changes of the 
crystal fields cannot be sufficiently accurate to convey 
the physical picture. Therefore, we consider a different 
analysis of these fields. In this approach, we calculate the 
centers of gravity of the id^^-y^ and t2g bands which are 
expected to be accurate due to their substantial shifts by 
about 2 eV caused by the crystal field. The obtained val- 
ues for A£2;2_„2 and Asxz+yz together with the similarly 
calculated Ap^ can be directly substituted into the equa- 
tion II]). For the known AEtot this equation provides an 
estimate for the minimal energy As^z^-r^ 

^^3z2-r2 ~ AEtot + 2^eff " Aex2-y2 
-4i°^^r^^lz^-r^,l - ^^xz+yz - Ap%, (5) 

which is sufficient for the stabilization of the orbitally 
reconstructed extended Cu 3d^ configuration. Here 
^3z2'Lr2 — —1-86 eV is the LDA-energy level of the unre- 
laxed 3(^322 _r2 orbital. 

In expression (O, the energies of the 3dx2^y2 bands 
can be represented as 

^x^-y^,a = ^''S-y'^.a + ^e// " <2_j,2_^^ . (6) 

In distinction to the reconstructed 3d3z2„^2 orbitals, the 
occupation numbers ri^2_y2 ^ of 3da;2_y2 states are af- 
fected only slightly by the decrease of h. In the 
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Ae^2°^^2 and the concentrations n'^2_y2 „ can be found 
from the centers of gravity of 3da;2„j,2 bands and by inte- 
gration of the densities of states obtained for the relaxed 
and optimized interfaces. Similarly to Ae3.2_j,2, the shift 
^£xz+yz can be easily determined from the changes of 
the LDA potentials e''*-"-' 
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FIG. 4: Band structure for different h where the top panel 
corresponds toh = 1.83 A and the bottom panel toh = 2.2 A. 
The size of the circles is proportional to the Cu 3(ij.2_y2 | 
character of the respective states. Here the coordinates of 
the symmetric points of the Brillouin zone are: A=(— vr, tt), 
R=(-7r, 0), and Z=(0,0). 



FIG. 5: Electron density map (in e/A^) in the energy range 
—4.5 eV < E < —2.5 eV below the Fermi level for unrelaxed 
and optimized interface distances. Here the electron density 
is plotted in the (x, z) plane with the a;-axis along the 0(1)- 
Cu-O(l) diagonal of the Cu04 plaquettes shown in the right 
panel of Fig. [1] 



latter property leads to small changes by 0.07 eV in the 
second local Coulomb contribution due to atomic relax- 
ation. Consequently, the effect of the interface relaxation 
on 3(1^2 _y2 is reduced to a modification of the LDA po- 
tentials 



= Ae 



(10) 



Ae 



(10) 



(7) 



where Ae'-l^'' 2 

Using the results of the LDA+U calculations, the shift 



1(0) 1 



x^ — y^^t7 x'^ — y'^^a' 



The band gravity centers calculated for Cu 3^2:2-^2, 
'idxzj^yz and for 2pz orbitals of 0(2) in superlattices with 
unrelaxed and optimized h are presented in Table[ll From 
this table, one can detect that a central feature of the in- 
terface relaxation is a strong shift of the center of the 
minority band 3c?a;2_j,2 down to the energy —1.07 eV 
below the Fermi level. The reason for this shift can be 
understood from the analysis of the band structures of 
the optimized and unrelaxed interfaces plotted in Fig. |4l 
In Fig. m the size of the circles is proportional to the Cu 
'id^2_y2^^ character of the respective states. The compar- 
ison shows a substantial Cu 3dj;2_j^2 ^ contribution to the 
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FIG. 6: Projected Cu 3d3^2_^2^|, 0(2) 2p^,| and Ti 
3d-j2^_j.2^l orbital densities of states at the interface between 
Cu02 and SrTiOs in the relaxed heterostructure (LSDA+f/ 
studies). The zero of energy is at the Fermi level. 



2p states of oxygens in the Cu02 planes which are located 
below the Fermi level. This contribution is especially sig- 
nificant at the energies about —(0.6 — 0.8) eV, at —2.5 eV 
as well as at about —(3.8 — 4) eV. Specifically, the inte- 
gration of the partial densities of states below the Fermi 
level shows an increase of the population An^2_y2^^ of 
Cu 3(ia;2_j,2 I levels by about 0.2 electrons in the relaxed 
structure. To check the numerical accuracy of this result 
we have also calculated the increase of the electron con- 
centration from equation ^ as An2.2_y2 ^ = n'^2_y2 ^ — 

n].2_y2 i = (£i2_j^2 1 — e°2_j^2 J - Ae^2L^2 |)/f/e// = 0.25 

with A(r^2'^],2 I = e'^2ly2^^ - e'^2ly2^^, which agrees weh 
with the value obtained by the direct integration. Be- 
sides the 3(ia.2_j^2 J state, the relaxation also results in an 
increase of the population of 'id,j.2_y2 '^ orbital by about 
0.1 electrons, as compared to the system with h — 2.2 A. 



Using the values of the band centers given in Table [H 
we find from Eq. ([5]) that the minimal energy shift 
Ae3°2_^2 sufficient for orbital reconstruction is given by 
Ae3^2_r2 ~ 1 eV. As the typical values of the crys- 
tal fields expected in this type of systems can approach 
2 — 3 eV, we can conclude that the mechanism of the 
orbital reconstruction realized through the stabilization 
of the extended Cu id^ configuration shown in Fig. [2] 
(bottom panel) appears to be a plausible scenario. We 
note that this reconstruction occurs despite a strong local 
Coulomb repulsion U — 8 eV which is overcome by the 
modifications of the structure of the Cu 3dx2-y2^^ and 
Sd3z2-r2,l bands. 



III. ELECTRONIC DENSITY AT ORBITAL 
RECONSTRUCTION 

To understand better the nature of the reconstruction 
of Cu 3(ij:2_j,2 |, we also analyze the (a;, z) electron den- 
sity maps (represented in Fig. [5]) in the energy range of 
the bonding states Cu 3dx2_y2-0 2p. Here the z-axis 
is directed perpendicular to (001). The decrease of the 
interface distance h leads to stronger electron repulsion 
in the nearest Cu02 and SrO layers. The effect of this 
repulsion is clearly visible in the bottom panel of Fig. O 
(h = 1.83 A) which shows electron depletion in the spa- 
tial region between Cu (located at a; = 2 A, z = 4.5 A) 
and oxygen 0(2) {z — 2.5 A) in the SrO-plane. The de- 
crease of the electron density along the z-direction has 
an extended character and involves Cu 3(i322_r2 |, 0(2) 
2pz^l and Ti 3d3z2_^2 ^ orbitals which may be interpreted 
as reduced hybridization between these orbitals. In Fig.[S] 
this decrease is manifested through a peak centered at 
E — 1 eV above the Fermi level which appears due to the 
splitting of the Sd^z^-r^ bands. The peak remains signif- 
icant even in the density of Ti 3(i3j,2_^2 ^ states which are 
more distant from the interface. The consequent increase 
of the difference between the electron occupancy of the 
majority and minority spin states also leads to the en- 
hancement of the magnetic moment of Cu from 0.65 fiB 
to 0.86 fiB^ an effect discussed in the previous section. 

In contrast to the electron depletion in z direction, 
in the Cu02-plane the interface relaxation results in a 
strong charge hybridization between Cu and 0(1) which 
is presented in more detail in Fig. [T] To estimate the 
modification of electron charge resulting from the com- 
bination of (i) electron depletion across the interface and 
(ii) electron accumulation in Cu02 planes, we have cal- 
culated planar (cc, y) and orthogonal (z) contributions 
to the electron density. To obtain the orthogonal part 
(i) (electron depletion), we have extracted the charge 
density in the range between 0.5 eV and 1.4 eV which 
corresponds to the new DOS-peak. The direct integra- 
tion of the obtained spatial density profiles results in the 
value rig ~ 1.45 carriers per unit cell. This value should 
be corrected by An^ on account of an overlap with the 
states 5dx2-y2, 2px and 2py which are shifted towards 
lower energies in the relaxed system (see Fig. [3]) . The in- 
tegration of the corresponding projected DOS parts gives 
Ang = 0.28. It should be noted that the latter integra- 
tion has been performed without consideration of extra 
charge in the interstitial regions which should give addi- 
tional contributions to An^. In our analysis, the inter- 
stitial contributions are expected to be small and the re- 
sulting Arig does not exceed the maximal value of about 
0.45-0.5. As a consequence, the obtained modification of 
the depleted charge density Anuf =n\ — Anf k, 1. This 
implies that approximately one electron per unit cell is 
depleted from the 3d3z2_r2 and 2pz orbitals of the in- 
terfacial bonding complex Cu-0(2)-Ti. In the interfacial 
stack Cu02-SrO-Ti02, about 70% of the depleted charge 
is located in the Cu02 layers which suggests a strong spa- 
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FIG. 7: Electron density map (in e/A^) in the energy range 
—4.5 eV < E < —2.5 eV below the Fermi level in (a;, y) 
plane for unrelaxed and optimized interface distances. Here 
z/c = 0.7. 

tial charge confinement to the 8(^322 orbitals of Cu. 

To calculate the amount of the electrons accumulated 
in the Cu02 planes of the relaxed structure (contribu- 
tion (ii)), we performed the integration of the hole charge 
from the Fermi level to the top of the 0(1) valence band 
which gives n^^ = 1 holes. These holes are located in 
the 2px.y states of 0(1) hybridized with the 3dx2-y2 or- 
bitals of Cu. The comparison with the value n^^{hi) — 2 
for the unrelaxed system provides the resultant amount 
rig" = 1 for the additional electron charge accumulated 
due to structural relaxation. It is remarkable that in their 
combination, both contributions (i) and (ii) give the total 
positive charge +2e required for the overall neutrality of 
the relaxed system. 



Therefore, in this mechanism for the orbital recon- 
struction one can distinguish two key factors which are 
responsible for the stabilization of the extended Cu 3(P 
configuration near structurally relaxed Cu02/SrTi03 in- 
terfaces: (i) the depiction of the electrons from a Sd^z'^-r^ 
orbital of Cu hybridized with 2p 0(2) and 8(^322 states 
of Ti, a process caused by the splitting of the 8(^322 _^2 
bands; (ii) the increase of electron charge density in the 
Cu 8(ia;2_j,2 orbitals which are strongly hybridized with 
the 2px^y bands of the oxygens 0(1) in the Cu02 planes. 
This increase is accompanied by a shift of the 3dx2_y2^^ 
band towards lower energies and it provides the energy 
gain sufficient for the orbital reconstruction. It is worth 
noting that the second stage (increase of electron con- 
centration on 0(1) 2p and Cu 3d^2_y2^^) in fact counter- 
balances the electron depopulation in stage (i). In this 
way, it minimizes the lost of energy due to the splitting of 
Cu 8c?322_^2 bands with the orbital reconstruction. The 
obtained charge self-regulation is in agreement with the 
negative feedback mechanism discussed in Ref. [l^ which 
supports its general character for the transition metal 
oxides. 



IV. CONCLUSIONS 

To summarize, we note that in the orbital reconstruc- 
tion mechanism, the effective change of the valency of Cu 
due to the splitting of 8^322 _^2 bands appears to be a 
plausible and realistic property. For the cuprate-titanate 
interfaces, the resulting electronic state is characterized 
by the electronic configuration 3d^^'^^ of Cu 3d orbitals 
extended by the changes of orbital occupancies of ions in 
the neighbouring SrO and Ti02 planes. Surprisingly, the 
strong hybridization between O 2px^y and Cu 8c?a;2_j,2 | 
states of Cu02 planes plays an essential role, making even 
moderate crystal fields acting on Cu 8^32-2 _^2 orbitals 
sufficient for the orbital reconstruction at Cu02-SrO in- 
terfaces and the consequential splitting of the 8^322 .^2 
states through the onsite Coulomb interaction U. An 
increase of the magnetic moment of Cu (from 0.65 
to 0.86 fiB) is predicted for the reconstructed interface 
state. 
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